Dynamics of magnetization coupled to a thermal bath of elastic modes 
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We study the dynamics of magnetization coupled to a thermal bath of elastic modes using a system 
plus reservoir approach with realistic magnetoelastic coupling. After integrating out the elastic 
modes we obtain a self-contained equation for the dynamics of the magnetization. We find explicit 
expressions for the memory friction kernel and hence, via the Fluctuation-Dissipation Theorem, 
for the spectral density of the magnetization thermal fluctuations. For magnetic samples in which 
the single domain approximation is valid, we derive an equation for the dynamics of the uniform 
mode. Finally we apply this equation to study the dynamics of the uniform magnetization mode 
in insulating ferromagnetic thin films. As experimental consequences we find that the fluctuation 
correlation time is of the order of the ratio between the film thickness, h, and the speed of sound in 
the magnet and that the line- width of the ferromagnetic resonance peak should scale as Bfh where 
B\ is the magnetoelastic coupling constant. 



I. INTRODUCTION 

Thermally induced fluctuations of the magnetization 
are responsible for one fundamental limit on the signal 
to noise ratio of small magnetoresistive sensors. 1 . The 
noise scales inversely with the volume of the sensors and 
peaks at frequencies 2,3 that are now close to the ever 
increasing data rate of magnetic storage devices. The 
increase of data rates combined with the continuing de- 
crease of the dimensions of the sensors makes magnetic 
noise inevitable and motivates work aimed at achieving 
a detailed understanding of its character. 

The standard approach toward modeling of magneti- 
zation fluctuations is to start from the Landau-Lifshitz- 
Gilbert-Brown equation 4 



dt M s 



5E u 



on 

an x — , 
at 



(1) 



where 7 is the gyromagnetic ratio, ft = M/M s is the 
magnetization direction, M is the magnetization, Ms the 
magnitude of the saturation magnetization, E the free 
energy and h a random magnetic field. This equation 
assumes that the characteristic time scale of the magne- 
tization dynamics is longer than the typical time scale 
of the environment that is responsible for the dissipative 
term proportional to a. In practice the use of this equa- 
tion is partially inconsistent, resulting in some practical 
limitations to its application^. The source of the prob- 
lem is that the dissipation is local in time. Because of the 
fluctuations dissipation theorem, this implicitly requires 
the random field to have white noise properties i.e. to 
have zero autocorrelation time. Since the contribution of 
the random field to the magnetization dynamics 7ft x h 
depends on ft, equation (Tj]) exhibits white multiplicative 
noise 7 . It follows that in order to integrate equationlflj 
reliably we need to track the evolution of ft on very short 
time scales for which the white noise approximation for 
h is likely to be unphysical. 

In this paper we address the physics that determines 



the correlation time of the random field. We start in sec- 
tion II by considering a formal model of a magnetic sys- 
tem coupled to an environment and specialize in Section 
III to an environment consisting of elastic modes. In Sec- 
tion IV we consider the case in which a single magnetic 
mode corresponding to coherent evolution of the magne- 
tization in a small single-domain system is coupled to the 
elastic environment. In Section V we consider a thin film 
geometry in which the magnetization is coupled to elastic 
modes of the system and its substrate. Finally in Section 
VI we conclude with a discussion of the possible role of 
other sources of dissipation, in particular dissipation due 
to particle-hole excitations in the case of metallic ferro- 
magnets. 



II. GENERIC RESERVOIR 

Calling q n the degrees of freedom of the reservoir, we 
consider the following form for the total Lagrangian: 

Qn 1 Qn 

]+£i[Cl(x),q n ]— A£[ft(x)], 
(2) 

where £s[ft(x), ft(x)] is the Lagrangian that describes 
the dynamics of the magnetization when not coupled to 
external degrees of freedom, Cn[q n , q n ] is the Lagrangian 
for the reservoir and £/[ft(x), q n ] is the interaction La- 
grangian that couples the magnetization to the reservoir 
degrees of freedom. The term A£[ft(x)] is a counter 
term that depends on ft and the parameters of the reser- 
voir but not on the dynamic variables of the reservoir^. 
This term is introduced to compensate a renormalization 
of the energy of the system caused by its coupling to the 
reservoir—. 

The Landau-Liftshitz equations for the decoupled sys- 
tem magnetization follow from the magnetic Lagrangian, 



Co = 



M 

—A [ft] -ft-£ s [ft] 



dx, 



(3) 



2 



where A is a vector field defined by the equation: Vjj x 
A[fi] = ft and 2?g[fi] is the magnetic free energy func- 
tional and Vm the volume of the ferromagnet. We model 
the reservoir as a set of classical degrees of freedom: 



(4) 



The Euler-Lagrange equations for the total Lagrangian 
J2J) yield the following coupled dynamical equations: 

d 

m n q n =- — [£ R (q n ,q n ) + Ci[Sl,q n ]] (5) 

^§ =nx irm [Es[n ' ™ ] ~ Cl[n -' qn] + AC[n]] - 

(6) 

When Ci is linear in the coordinates of the bath, we can 
formally integrate (J5J to get q^ n \t) as a function only 
of the initial conditions and O and then insert the re- 
sult in jnj to eliminate the reservoir coordinates from 
the dynamical equations for H, integrating out the reser- 
voir degrees of freedom. An example of the application 
of this procedure for a quantum mechanical model of the 
interaction between magnetization and reservoir degrees 
of freedom can be found in Ref. ^3 



III. MAGNETIZATION COUPLED TO ELASTIC 
MODES: GENERAL 

If we consider only long wavelength vibrations we can 
treat the lattice as a continuous medium and use results 
from elasticity theory. The potential energy functional, 
Er, of the elastic medium can then be expressed in terms 
of the strain tensor Uj 



1 / diM 

2 [dx, 



dxi 



where u is the displacement vector field. 

We want to study the dynamics of the magnetization 
when coupled to elastic deformations of the systemii. 
We will be interested in applying our results to polycrys- 
talline elastic media which can be treated as isotropic to 
a good approximation. (It's quite straightforward, albeit 
quite tedious, to extend our results to the case of non- 
isotropic media with specific lattice symmetries). For 
isotropic elastic media it follows from general symmetry 
considerations that, to lowest order, we can express the 
magnetoelastic energy in the form^ 



E I = B l Y J 



QifljUij dx 



(7) 



where B\ is the magnetoelastic coupling constant. For 
the case of soft ferromagnet thin films, the main contri- 
bution to the magnetoelastic energy will be given by the 



magnetostatic energy dependence on strain. This contri- 
bution to Ei is normally referred as the the form effect. 
The constant B\ can be extracted from magnetostriction 
data. For an isotropic elastic medium with isotropic mag- 
netostriction, A, we havei£ that 



* = Ix, E 



2 2-cr' 



(8) 



where E is the Young's modulus and a the Poisson's 
ratio. 

The Lagrangian for an elastic reservoir Cr is, 



Cll = \ J P™ 2 ^- E R-- 



(9) 



where p is the mass density, V the total volume of the 
elastic medium (magnetic film plus substrate) and Er is 
given bjiii 



Ee 



E x 2 lE x ^ 

2(1 + a) .4- Ul1 + 2(l + a)(l-2a) 2-r ' 

V ' % J — 1 V ' K ' 1=1 



dx 



The equation of motion for the displacement will then 
be, 



(id 



It will prove useful to expand u in terms of the elastic 
normal modes f ( n ' : 



u = ^g (n) (t)f (n) (x) 



(12) 



where the functions f satisfy the boundary conditions 
appropriate for u and satisfy: 



<pfW(x); fieN 
S nrr 



5f(")(x) 

± J pf(")(x)-f( ro )(x)dx 



(13) 
(14) 



where M is the total mass, M = J v pdx. 

In terms of the degrees of freedom, q^ n \ we have: 

d = -Ej = -B 1 J2 q (n) E / n A'4 n) dx ( 15 ) 



with 



P (n) 



I (dft ] 9f, 



dxi 



dxi 



We then see that the interaction Lagrangian is linear in 
the coordinates c^™- 1 , with coupling constants: 



c W[n] = W <),<),/,;'dx 

i,3 JV ™ 



(16) 
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This property will allow us to integrate out the reservoir 
degrees of freedom to obtain an equation for the dynam- 
ics of the magnetization in term of ft alone. 

Let's first discuss the dynamics of the reservoir degrees 
of freedom gW, Using equations l)ll [l -l)14 [l we find the 
dynamical equations: 



(17) 



Integrating l|17(l we find 



n( n ) I 

g W(t) =g(«)| cos(a; n t) + s in(u;„t) 



_Bi_ 

Mv r , 



sm(oj n (t-t'))c^[fl(t')]dt', (18) 



of the magnetization we find: 
on t SE S 7 SACjfl) 



M, ^ 



Bi 



M s Sfl 



Sfi 



x ( 9 (n) | cos(w„<) + 



g (n) |o 



sin(w ra t) 



c W[n(t)] 



^ C («)[n(0)]cos(^) w 



x,/ 



Bl ^ dt' 



Sc^ 

dx' cos(w„(t — i')) 



sn 



dfl 







x , SSI 


x,i 



(21) 



The counter term AC of the total Lagrangian is defined 
to cancel the frequency renormalizing term: 



where q^lo and q^lo are the initial values of and 
q(n) respectively. The coupling of the magnetization to 
the reservoir will cause damping and frequency renormal- 
ization. In order to be able to separate the two effects is 
useful to integrate the last term on the right hand side 
of i(T%|) by parts obtaining: 



Bl 

Mufi 



[«(*)] 



Sc^ 



Sfl 



It follows from Eq. JTBI) that 

Set' 



sn. 



Q f (n) g f (n) 



dx,. 



dx. 



(22) 



(23) 



/}(«) I 

g («)(t) =g(")| cos(w„t) + 2—!° sinKt) 



dt' 



u n JO 

Sc^ 



V M 



Sfl 



cos(uj n (t — t')) 
dfl , 



(19) 



To simplify and extract the physical content from these 
cumbersome equations, we identify the memory friction 
kernel tensor 7j- m : 



ijm 



(M',x,x')^e(t-i')E 



M s Mui 
ScS 



2 cos(w„(i - 1')) 



Sfl, 



ScW 



c't' SQ 3 



x./ 



(24) 



where 0(i — t') is the Heaviside function. We also recog- 
nize the random field h: 



Using the expression of the interaction Lagrangian given 
by (jT5|) and the definition of the coupling constants 
we have: 



Sfl ^ q Sfl 



(20) 



Combining equations © , (|19fl and l|2U|) for the dynamics 



M, 



{ 1 gC") [ ] <5 c ( n ) 

9 w |ocos(w„t) H sin(w n t) 

Assuming that the distribution of initial positions of the 
environment degrees of freedom follows the canonical 
classical equilibrium density for the unperturbed reser- 
voir wc find that 



(hj(x,t)h m (x',t')) 



h(x,i)> = 0, 

2K B T 



1 M S 



(26) 

7jm (M',x,x'). (27) 



4 



In terms of 7 jm and h the dynamical equation for O takes 
the form: 



d^li „ 7 SE S 



rt 



eijifli / dt' / dx'y^jj m (t,t',x,x 
Jo Jv M 



, dO r , 



dt' 



7 B l - -^ c(Tl) 



i.n 



50 < 



The final term is an artifact of the assumption that in 
the initial state the reservoir was decoupled from the 
system^ 5 . Dropping this term, the dynamical equations 
for magnetization coupled to a thermal bath of elastic 
modes is: 



dO t 7 5E S 



rt 



CijiOi / dt' / dx' y^7-, m (t,t , ,x, x 
Jo Jv M 



, dOj, 



dt' 



(28) 



with 7j m defined by (|24() and h a random field with sta- 
tistical properties given by |23) and J2ZJ|. Equation QgJ 
is quite general. In particular notice that to obtain Ij28(l 
we didn't perform any expansion in O. As a consequence, 
as long as we keep the exact form for Eg(£l), equations 
(I28JI includes also the effects of spin wave interactions. 
In principle we could also include in Es a term to take 
into account the scattering of spin waves due to disor- 
der. Equation (|28(l does not, however, take into account 
the coupling between the magnetization and particle-hole 
excitations. As we discuss in Section VI, this coupling 
appears to be of critical importance in many metallic 
ferromagnets. 

Equation (|28(l is very different from the standard 
stochastic Landau-Lifshitz-Gilbert (s-LLG) equation, 
Eq. QJ. Because the magnetoelastic energy, E], J7|), is 
nonlinear in the magnetization, in l|28|) both the damping 
kernel and the random field depend on the magnetization 
and therefore are state dependent. This is in contrast 
with the s-LLG equation for which both the damping 
kernel, a5(t — if), and the random field are independent 

of n. 

Another difference between Eq. (|28|l and the s- 
LLG equation is that the damping kernel, "fj m , is in 
general a tensor. The tensor character of the damp- 
ing has been suggested previously on phenomenological 
grounds^. Starting from the physical coupling Q, in 
our approach the tensor character of 7 jm appears natu- 
rally as a consequence of: (a) the nonlincarity in O of 
the magnetoelastic coupling (JJJ, (b) the anisotropy of 
the elastic modes due to the boundary conditions and/or 
anisotropy of the elastic properties. For small oscillations 
of O around its equilibrium (up to quadratic order), the 
kernel 7 jm can be assumed to be independent of O. Even 
in this linearized case, the damping kernel that appears 



in H28fl will still have a tensor form due to the anisotropy 
of the elastic modes. 

As mentioned above, the standard s-LLG damping ker- 
nel is simply a5(t — t'), i.e. the damping is frequency- 
independent. As a consequence, from the Fluctuation 
Dissipation Theorem, we have that the spectrum of 
the random field that appears in Q is also frequency- 
independent. This differs from equation (|28|l for which 
the damping kernel, and therefore the spectrum of the 
random field, is frequency-dependent. 

Given the geometry and the material properties of the 
system we can find the elastic modes, f( n >, and then in- 
tegrate equation l|28l) using a micromagnetic approach. 
The integration of equation l|28|l could give insight in 
particular on the damping of the uniform magnetization 
mode for different geometries and show the range of va- 
lidity of the classic picture^ of a two stage damping pro- 
cess in which the motion of the coherent magnetization 
induces non uniform magnetic modes on short time scales 
that then decay to lattice vibrations. 



IV. MAGNETIZATION COUPLED TO ELASTIC 
MODES: UNIFORM MAGNETIZATION 

We now study the dynamics of the uniform magnetic 
mode in the case when we can neglect its interaction with 
spin waves and the only coupling to external degrees of 
freedom is magnetoelastic. Projecting Eq. I|28|) on the 
uniform mode we find that 



f SE S 7 f 

fi,-L [dt> /dx fdx'^ V (M',x,x')^ 
Vm Jo Jv m Jv m m at 

(29) 



Let's define the space averaged error field 



1 



h(i) = — / h(x,t)dx, 

VM Jv M 



the damping kernel 

Ijm (*,*') = 77- / dx / dx'7 jm (£,i',x,x') 
Vm Jv m Jv m 



and the coefficients 



» 



Vm 



8cW 
~STh 



dx. 



Using the fact that f2 is uniform we obtain 



» 



(n) 



dxi 



dxi 



dx. 



(30) 



In terms of the coefficients c, we can then write: 



hi = 



B 



M S V M 



.(«) 



i ^ q (n) \o N 

q {n > | cos(uj n t) H sin(u; n t) 
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and 



j jm = G(t-t') 



M S MV M ^ u>. 



The uniform magnetization dynamics can then be ex- 
pressed in terms of the spatially averaged random field h 
and memory friction kernel jji : 



d^ _ f SEs 

At eijl l V M M s J Vm SQj 



dx + jtijiClihj 



with 



and 



{hj{t)h m {t')) 



2K B T 
jV M M a 



Jjm(t,t'). 



(32) 



(33) 



(34) 



V. THIN FILM UNIFORM MAGNETIZATION 
DYNAMICS 

We now apply equation l|32|l to study the dynamics of 
the uniform magnetization in a thin ferromagnetic film 
placed on top of a non magnetic substrate and covered 
by a non magnetic capping layer, as illustrated in Fig. ^ 
We assume that all media are polycrystalline and treat 
them as isotropic. We will assume the lateral size, L s , 
Fig. ^ to be much bigger than the film thickness h. 
Notice that if we take L s bigger than the domain wall 
width our assumption that the non uniform magnetic 
modes are quenched wouldn't be valid anymore. We will 
consider only oscillations of the magnetization around 
an equilibrium position parallel to the X3 axis so that we 
can calculate the damping kernel tensor "fj m assuming 
the elastic modes to depend only on £3. Otherwise, to 
find the correct damping kernel, we would have to take 
into account the fact that the lateral size, L s , is finite 
and solve the full 3D elasticity problem for the elastic 
modes. 



A. Damping kernel and random field 

To find the dynamics of the magnetization using equa- 
tion l-'il'l) we need to evaluate the memory friction kernel 
7j m . The first step in this calculation is the determina- 
tion of the elastic normal modes f( n ) which satisfy the 
following equation: 



(n) 



E 



2(1 



_y2 f (n)_ 



E 



2(l + o-)(l - 2a) 



V(V-f (n) ). 



Substrate (S) 



Y/ ■ 

Y . 

"1 x. 



FIG. 1: Geometry considered for the case of a thin ferromag- 
netic film on a non-magnetic substrate. 



We allow the film, the substrate, and the capping layer to 
have different elastic properties and solve equation 135|) 
separately in the different subsystems using the appro- 
priate elastic constants. We assume for the sake of def- 
initeness that the substrate and capping layer material 
is identical. We then match solutions by imposing the 
continuity of displacement and stresses at the interfaces 
X3 = 0, and X3 = h. As boundary conditions we assume 
the top surface of the capping layer to be free and no 
displacement at the bottom of the substrate. 

Because in our case the elastic modes only depend on 
X3, Eq. l|3Tjjl simplifies to 



» 



1O3H 



with 



(n) 



ft\h)-fr(o). 



in). 



The spatially averaged damping coefficients have a simple 
expression in terms of the A/^ : 



T 2 r2 

,J 1 ; Mh 



[Af(»)]2 

J2 1 \ J cosK(t-O) 



^7, 



X [<%0 3 (*) + fii(t)*3i]fci«3(f) + fii(i%i]. (36) 

Eq. 1201 follows from the completeness relation of the po- 
larization vectors. Once we know the coefficients Af( n \ 
Eqs. I|36|l . (|33|l . and Ij34(l completely specify the dynam- 
ical equation l|32(l for the magnetization. 

As an example we consider the case of a polycrystalline 
ferromagnetic thin film, like YIG, placed on a substrate 
of a polycrystalline paramagnet like Tantalum, Ta. As 
typical values we takei£ the ones listed in Table [I] For 
the magnetostriction we assume A = 2 x 10~ 6 . Using 
equation JSJ, we find that B\ = 4 x 10 6 ergs/cm 3 . Given 
the elastic modes implied by these parameter values, we 
can calculate the coefficients Af^ n \ Once we know the 
coefficients A/-*-"-* we have all the elements to completely 



(35) specify equation (|32p. 
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M s Mh 



TABLE I: Elastic properties, ct, ci are the transverse and 
longitudinal speed of sound respectively. 





Magnetic Film 


Substrate/Capping Layer 


E 


200 Gpa 


180 Gpa 


a 


0.33 


0.33 


P 


5.0 g/cm 3 


16.6 g/cm 3 


Ct 


4.0 km/s 


2.0 km/s 


Cl 


5.0 km/s 


4.1 km/s 



We generate a stochastic field h with the correct sta- 
tistical properties by using its Fourier representation. To 
obtain 



we chooseASi 



where 



(y(t)y(t')) = G(t-t>) 



i r°° 

I/H = 5-/ yit)e-^dt 
*n J-oo 



(37) 



(38) 



and 



G{w) = 



1 

2V 



G(r)e" 



dr. 



In our case we have from Eq. I|31[) . that the memory 
friction kernel 7,7 depends separately on t and t' . As 
aconsequence, through lp"""*) . we have that the average 
(h(t)h(t')) doesn't depend only on the time difference 
t = t — t'. The random field h(t) therefore doesn't de- 
fine an ergodic process and in particular we cannot use 
equation (|38|) . For this reason it is convenient to define 
the auxiliary random variables: 



i \ q {n) \o 

q (n) o cos(w„i) H sin(a; n t) 



and the auxiliary kernels: 

[A f(™)l2 

9i (t - 1') = e(t - o V [ Jt 2 J cos( w „(t - 0) 



so that we have: 

(x i (t)x j (t')) = ?- 1 ^-g i (t-t')6 ij . 

The random variables Xi (t) therefore describe an ergodic 
process and we can use equation (13*51 to generate them. 
In terms of Xi and gi we have: 



h, = 



D 



MM 



5^5i[<5tfn 3 (t) + ni(t)« 



33 J 



(39) 



To generate the random field and calculate 7^ we then 
have to calculate the quan tities ffj(r ) and their Fourier 
transforms gi(to). Figures 2(a)|3(bj show some typical 
profiles for g%{r) and gi(w) using for the mechanical prop- 
erties the values of table"]] We find that in general gi{r) 
doesn't depend on the thickness of the capping layer V . 

In the limit in which we can linearize the magnetoelas- 
tic interaction with respect to O, we have: 



M.Mh 



9 3 ( T ) S . 



jm- 



(40) 



The damping kernel is diagonal with components equal, 
apart from an overall constant, to gj{r), in contrast to 
the s-LLG equation for which we have 7,7 (t) = a5{r)8ji. 
The power spectrum of the random field component, 
hj, is then proportional to gj(to), in contrast to the s- 
LLG equation for which the power spectrum of each 
component hj is simply a constant. Notice that even 
in this limit jj m preserves its tensor form due to the 
anisotropy of the elastic modes. In our specific case we 
have g\ — g2 ^ 53 due to the difference between the 
transverse and longitudinal speeds of sound. 

From figures 2(a) 2(b)| we see that gi(j) goes to zero 
for times longer than T£> « 5 10~ 2 t = hh/c, where c = 
c t .M is the transverse speed of sound in the magnet. For 
a film 20 nm thick we then find td « lOps . When the 
relevant frequencies of £1 are much lower than 1/tjj, we 
can replace the damping kernel given by l|40|l with the 
simple kernel 



Jjm = ljosS(T)6j 



with Tjeff given by: 



7jeff 



M*Mh 



(41) 



In this limit we recover a damping kernel of the same form 
as the one that appears in the s-LLG equation. Here 7 JO ff 
is the equivalent to a in J*Q). From the results shown in 
Fig 2(a) 2(b) we see that we have 



and then: 



Tjeff 



h 2 (L + h + L') 



lB\h 



M sP c 3 



(42) 



We find that the damping of magnetic modes in thin 
films is proportional to B\h. Assuming the values given 
in Table ITT1 we find 7i e g = 7 2o ff ~ 2 x 10~ 4 . 
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TABLE II: 



Quantity 


Value 


7 


1.76 x 10 7 s- 1 G- 1 


Bi 


4 x 10 6 ergs/cm 3 ; 


M s 


150G 


L 


lyxa 


ft 


20nm 




h = 0.01 L 
h = 0.02 L 



0.1 



(a) 




h = 0.01L 

h = 0.02L 



-0.1 



0.1 



(b) 



FIG. 2: (Color online) Profiles of <?i = gi (r)c 2 /[h(L + h + L')}, 
(a), and §3 = (?3(r)c 2 /[/i(L + ft + £')]> C 5 )) f° r the case °f a 
thin magnetic film on a Tantalum substrate; to = L/ct,M- For 
the standard s-LLG equation (r) would simply be a Dirac's 
delta centered at r = 0. 



0.05 



0.04 



,3 

<^§' 0.03 

a 
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i 0.02 - \ 



0.01 - 



• h = 0.01 L 

• h = 0.02 L 



V 




600 



800 



(a) 



0.015 



0.01 - 



• h = 0.01 L 

• h = 0.02 L 



3 



0.005 




(b) 



FIG. 3: (Color online) Values of Re[&(u>)] = 
Re[ ff i(w)]c 2 /[ft(L + ft + I/)] (a) and Re[g 3 (tu)] = 
Re[g 3 (io)]c 2 /[h(L + ft + L')] (b) at the elastic modes 
frequencies {u>„} for the case of a thin magnetic film on a 
Tantalum substrate. Shown are the values for ft = 0.01L, 
diamonds, and ft = 0.02L, circles. For any ui„ Re[gi(wn)] is 
unique even though this is not completely evident from the 
figure because in order to show the behavior of the auxiliary 
kernels over a wide frequency range, the resolution is not 
high enough to always show the separation between the 
single points. For the standard s-LLG equation gi(u>) would 
simply be a constant. 
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B. Integration 

After generating the random field h in the way de- 
scribed above we can proceed in integrating equation 
(J32J- We assume SE S /Sfl = -V M M s H cS with H cff = 
(0,0,-ffeff) and H e g simply a constant. Let's define the 
dimensionlcss quantities: 



t = jH cS t; H 



H 



off 



off 



H, 



cff 



H, 



cff 



Jjm = —rr~ ; 1 = 



7 ff ff' ~ H cS M s V M ' 
then in dimensionless form equation 132(1 takes the form, 
d^z 



df 



dflQi f di'J^Jjmdi')^ (43) 



with 



(fy(t)> = 0; (h 3 {t)h m (i')) = fj jm (t, t'). (44) 

Similarly, for 5Es/5Q = — VMM s H c g, the standard 
s-LLG equation, (I}, for the uniform mode, takes the 
dimensionlcss form: 

df2 d$~i 
— r = - ft x H eff + $7 x h + atl x — - (45) 
dt di 



with 



{h 3 (t)) = 0; (hj(i)h m (i')) = afS{t - t'). 



(46) 



the expression l)39|l and for gi(r), <7i(o>) 
2(a)|3(b)| and assuming T 



Using for j jr 
the results shown in figures 

1CP 2 and the values given in Table IIII we integrate equa- 
tion (|43|l . We used the stochastic Heun scheme that 
ensures convergence to the Stratonovich solution even 
in the limit of zero autocorrelation time for the ran- 
dom field 7 . The results of the integration are shown 
in figures 4(a)|4(b)|5(a) As initial condition we took 
ft = (0.6, 0,0.8), dS 2/ dt = 0. 

We then integrated equation ((45(1 setting a = 7i e g with 
7icff calculated using 141(1 . The results of the integration 



are shown in fi gures |4(a)|4(b)|5(b)| 

From figures |4(a)|5(b)| we see that on average equation 
(|43(l and ((45(1 give very similar results. This is expected 
because for the initial conditions chosen we are in the 
limit of small oscillations around the equilibrium position 
and therefore the dependence of 7j m on Jl is negligible. 
The main differences, for the case considered, between 
the results obtained using 1(43(1 and ((45(1 are in the random 
fluctuations of ft. This is a consequence of the different 
correlation in time of the random field h used in l(43|) and 
((45(1 . For example we notice that equation ((43(1 seems to 
give a less noisy dynamics than ((45(1 even though for both 



simulation \h\ 2 is of the same order of magnitude. If we 
zoom on a short time interval, fig. |4(b)| as a matter of 
fact, we see that on very short time scales the amplitude 
of the random fluctuations for the two simulations is very 
similar. However for 145(1 fluctuations with the same sign 
are much more likely than for 1(43(1 . This is due to the 
different spectral density of the random field. For I45|) 
we simply have |/ij(o;)| 2 = aT, whereas for 1(43(1 |/jj(o;)| 2 
is equal to gj(u>) (considering that for our simulation, to 
a good approximation, we can neglect the dependence of 
the random field onfi). In particular for ((43(1 hj (uj) \ 2 has 
a low frequency cutoff at lu — u>q = c t ,M / ' L, where c^m is 
the transverse speed of sound in the magnet. This implies 
that for 1(43(1 we have a much lower probability than for 
((45(1 to have consecutive fluctuations of the random field 
with the same sign with the result that the dynamics 
appears less noisy. 



VI. DISCUSSION AND CONCLUSIONS 

In this paper we derived the equation for the dynamics 
of the magnetization taking into account its coupling to 
the lattice vibrations. The equation that we obtain, 1(28(1 . 
is quite general. Equation ((28(1 will have the same form 
also if we include spin-spin and spin-disorder interactions. 
To take into account these phenomena it is necessary only 
to add the appropriate terms to the energy functional 

E s [n}. 

From the general equation we derived the equation, 
1(32(1 . for the dynamics of the uniform magnetic mode in 
a thin magnetic film when nonuniform magnetic modes 
can be assumed frozen out. We find that in general 
the random field that appears in the dynamical equa- 
tion for the magnetization has a correlation time, to, 
of the order of the ratio between the film thickness, h, 
and the sound velocity c. When the timescale for the 
dynamics of the magnetization is much longer that To, 
we recover the stochastic LLG equation. In this limit 
we calculated the value of the effective Gilbert damping 
constant, a. For typical ferromagnetic insulators, like 
YIG, we find a « 10~ 4 , in good agreement with the val- 
ues measured in experimentsi&IS. We can then conclude 
that for magnetic insulators magnetoelastic coupling is 
the main source of magnetization damping. 

Our work predicts that magnetic resonance experi- 
ments on ferromagnetic insulators should be able to ob- 
serve the anisotropy of the damping and as a consequence 
of the correlation of the thermal fluctuations. With our 
theory is possible to exactly calculate the spectral den- 
sity of the thermal fluctuations. The spectral densities 
for small samples will be different from the one observed 
in bulk experiments because of the discreteness of the 
elastic modes. It would be very interesting to test these 
results with new experiments on small ferromagnetic in- 
sulators samples. In particular for thin films one exper- 
imental consequence of our work is that the correlation 
time of the magnetic fluctuations will be of the order of 
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cf 




0.78 1 1 1 1 1 1 1 

2000 A 4000 6000 

t 



(a) 



0.996 




993 I ' i ' i ' i ' i ' 

4685 4686 4687 A 4688 4689 4690 

t 

(b) 

FIG. 4: (Color online) ^3 as a function of time obtained 
integrating the standard s-LLG equation, (I45L and equation 
l|4i-{fl . In (a) the trace obtained using equation 1)43(1 has been 
offset up by +0.05 for clarity. In (b) the trace of Q3 is shown 
on a short time scale. 



h/c where h is the thickness of the ferromagnetic film 
and c the speed of sound in the magnet. We also found 
that in the limit when the magnetization evolves on time 
scales much bigger than h/c the damping of the magnetic 
modes is directly proportional to B\h. The line width of 
the ferromagnetic resonance peak in insulating ferromag- 
netic thin films should therefore scale as Bfh, which, in 
principle, can be confirmed experimentally. 

For ferromagnetic metals, like permalloy, we also find 
a w 10~ 4 . This value is about two orders of magnitude 




-1 1 1 1 • 1 • 1 

2000 A 4000 6000 

t 

(b) 

FIG. 5: Envelope curves of the trace of fh in time as obtained 
integrating equation 1431 (a) and equation 14511 (b). f2i os- 
cillates between the maximum and minimum value given by 
the envelope curves with frequency yH c g, equal to 1 in the 
dimensionless units used. 



smaller than the value observed experimentally 20 . The 
reason is that in ferromagnetic metals the electronic de- 
grees of freedom are the main source of dissipation for the 
magnetizationSi* 2 ^. Starting from a model of localized d 
spins exchange-coupled to s-band electron, the interac- 
tion Lagrangian will be: 

Ci = J sd J dxfi(x) • s(x) 

where J sc i is the exchange coupling constant and s is the 
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conduction electrons spin density: 

S W = ^H*i( x Kfc* b (x) 



where $ are the s-band carrier field operators and T ab the 
representation of the spin operator in terms of Pauli ma- 
trices. By integrating out the s-band degrees of freedom, 
in the linear response approximation Sinova et al^, for 
the damping of the uniform magnetic mode find: 



a = hm ' 



t^o 2M s Hu J (2tt) 



de 



^|(^(k)|r+|V^(k))| 5 



x / —A aM (e)A bM (e + fuu)[f(e)-f{e + huj)} (47) 

where A a ^(e) and Ab^(e) are the spectral functions for 
s-band quasiparticles and /(e) is the Fermi-Dirac distri- 
bution. Equation (|47|l gives zero damping unless there is 
a finite-measure Fermi surface area with spin degeneracy 
or there is a broadening of the spectral function due to 
disorder 24 . Characterizing the quasiparticle broadening 
by a simple number T = h/r s , where r s is the quasipar- 
ticle lifetime, we can assume: 



-4a,k(e) 



(e-£a,k) 2 +r 2 /4- 



(48) 



Inserting this expression for the spectral functions in 
B47J1 we find a as a function of the phenomenological 
scattering rate T. Notice that l|47|l includes the contri- 
bution both of intra-band, and inter-band^SSiSI quasi- 
particles scattering events. The intra-band contribu- 
tion is due to spin-flip scattering within a spin-split 
band and is nonzero only when intrinsic spin-orbit cou- 
pling is present. From equation (|47fl . using the ex- 
pression for A a: k(e) given in (|4*5|l . we see that in the 
limit of weak disorder, small T, the intra-band contri- 
bution to a is proportional to 1/r, in agreement with 
experimental results for clean ferromagnetic metals with 
strong spin-orbit coupling2£i2L2fli2i and previous theoret- 
ical wor k 2l ?i 2 7i? 2 i??i? 4 i?5 . Similarly from (gZj) we sec that 
the inter-band contribution to a is proportional to T. 
This result agrees with the experimental results for fer- 
romagnetic metals with strong disordej2& and previous 
theoretical workSSiSSiSl. Notice that equation 1)47(1 im- 
plicitly also includes the contribution due to the so called 
spin-pumping e&eciM£££2JiL2L in which spins are trans- 
ferred from the ferromagnetic film to adjacent normal 
metal layers as a consequence of the precession of the 
magnetization. In order to calculate this effect in first 
approximation we simply have to substitute in (|47|l the 
conduction band quasiparticle states, ip, calculated tak- 
ing into account the heterogeneity of the sample. Assum- 
ing for the scattering rate, l/r s , typical values estimated 
by transport experiments, from equation (|47|l we find val- 
ues of a in good agreement with experiments. 

In summary we have studied in detail the effect of the 
magnetoelastic coupling to the dynamics of the magne- 



tization. Starting from a realistic form for the magne- 
toelastic coupling we have found the expression for the 
damping kernel, Jj m . We find that in general jj m is 
a nondiagonal tensor non-local in time and space. The 
knowledge of the exact expression of jj m allows us to cor- 
rectly take into account the autocorrelation of the noise 
term overcoming the zero correlation approximation of 
the stochastic Landau-Lifshitz-Gilbert equation. We find 
that for thin films for which the single domain approxi- 
mation is valid, both the damping and the fluctuations 
correlation time are proportional to the film thickness. 
Our results apply to systems for which the direct cou- 
pling of the magnetization to the lattice vibrations is the 
main source of the magnetization relaxation. We have 
shown that this is the case for ferromagnetic insulators 
whereas for ferromagnetic metals the magnetization re- 
laxation is mainly due to the s-d exchange coupling. 
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APPENDIX A: SIMPLE ESTIMATE OF -y eS 

Let's start from the definition of jj m (equation (|31Jl ): 
l 3 m = e^-OAr^^cf ( t )cW(Ocos( Wn (t-f)) 

n « 

(Al) 

where A\ = r )B\jM s MVM- For the case of thin film we 
found: 



where: 



A/fWf n) W-/i tt) (0). 



An) 

Notice that, by definition, are dimensionless and so 

are the quantities A/} . Assuming that at equilibrium 
is fl = (0, 0, 1) and keeping only the leading terms in fi 
in the expression for c, we have: 

4 n) =L 2 A// n) . 

Let's now expand the collective index n in its compo- 
nents: k, s where s is the polarization index of the elastic 
modes. Then, using the completeness of the polarization 
vectors and the fact that the polarization directions are 



11 



parallel to the axis Xi,X2,X3 we have: 

J jm =&(t t')A x ^c { ;\t)ct ] (f) cosK(i - 0) 
=Q(t - O^L^ ^-A/^A/ r k / cos( Wk , s (i - 0) 
=Q(t - t^AiLiY-LAftAfaSjmCos^jit - 0) 

k W k,j 

Now note that: 

M = pL 2 s L(l + h + L'); V M = L 2 h: 

where h = h/L,L' = L 1 / L. Then we can write: 

l jm =Q{t-t')- 



where ojq — c/L. In this approximation we can then 
write: 

l 3m {T) « Q(t - t') -^3 1 -L^e— DT . 

M s pL(l + h + L')h u 2 D uj 

Integrating this expression between r = and r = oo we 
find 



7eff = 



IB* 



1 1 



M s pL(l + h + L')h u D 
1 B 2 h 2 L 



M s pL{\ + h + L')h c c 



7^i 2 h 



iBl 



M sP L{l + h + L')h 

x ~T~ S Jm cos(wk,j(* - *'))• 
k w kj 

For small enough /i/L we can assume A/j fe ^ « kh with a 
cutoff for krj such that kr>h — l.We can then define the 
cutoff frequency lud = cko = c/h. With this approxima- 
tion we have: 

2^ ~~T cos ( w k,j(£ - t )) 



M s p(l + h + L')c 3 ' 
To be more accurate let's define the functions 



(A2) 



1 c 2 ^ [A/jf 

..2 c °s(^k,i(r)) 



h(l + ft + L') L k W L 
so that we can write 

The functions <?j(t) are plotted in figure 2. Integrating 
<7j(r) between and oo we find: 



= J2 ,,2 \ ,,2 COSK,, {t - t')) 



UJ 



D JO 
1 1 



S(w - w k — T cos(uj(t - t')) duj anc j then finally 

UJ Z + UJ 



UJ 



D 



Ujfj UJ Jo + oj d 

1 1 

:— o woe 

0J D UJ 



jI+ujI 

-U) D (t-t') 



cos(u>(t — t')) duj 



7cff 



IB 2 h 



M s pc 2 c ' 

analogously to what we found previously, (JA2I 
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